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Abstract — A steady two dimensional mixed convective 
magnetic hydrodynamic flow over a porous vertical plate has 
been considered. The governing equations and boundary 
conditions are non dimensionalized then transformed by one 
point group transformation and finally solved numerically by 
using Runge-Kutta fourth-fifth order numerical method with 
shooting technique. Effects of radiation parameter, suction/ 
injection, porosity, magnetic parameter, Schmidt number and 
Prandtl number on velocity, temperature and concentration 
profile are discussed and presented graphically. It is observed 
that velocity and temperature decrease while concentration 
increases with an increase in radiation parameter. The values of 
physical quantities skin friction coefficient, Nusselt number and 
Sherwood number are also tabulated with the variation of 
physical parameters. 

Index Terms — Heat transfer, MHD, Mass transfer, porosity, 
Radiation 


I. Introduction 

The effects of radiation on a MHD fluid flow due to a 
stretching/ shrinking surface have great interest of researchers. 
Chamkha [2] investigated the coupled heat and mass transfer 
by natural convection of Newtonian fluids in the presence of 
magnetic field and radiation effects. Mahmoud [9] presented 
a study of the flow and heat transfer of an incompressible 
viscous electrically conducting fluid over a continuously 
moving vertical infinite plate with uniform suction and heat 
flux in the presence of radiation taking into account the effects 
of variable viscosity. Cortell [3] presented a numerical 
analysis for flow and heat transfer in a viscous fluid over a 
sheet nonlinearly stretched with effects of thermal radiation. 
Radiation effect on fluid flow with and without applying a 
magnetic field have been investigated among others: Bataller 
[1], Hsiao [6], Jat and Chaudhary [7] etc. Mahanti and 
Gaur [8] investigated the effects of viscosity and thermal 
conductivity which vary linearly on steady free convective 
flow of a viscous incompressible fluid along an isothermal 
vertical plate in the presence of heat sink. Seddeek et al. [10] 
have presented similarity representation of MHD flow with 
heat transfer taking into consideration variable viscosity and 
thermal conductivity. 
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Hamad et al. [5] studied a steady laminar 2-D MHD 
viscous incompressible flow over a permeable flat plate with 
thermal convective boundary condition, radiation effects and 
similarity representation of the governing partial differential 
equations obtained by group method. Ferdows et al. [4] 
investigated a free convective heat and mass transfer flow over 
a moving permeable flat vertical stretching sheet in the 
presence of non-uniform magnetic field. Effects of thermal 
radiation and convective surface boundary condition on steady 
boundary layer flow of a viscous incompressible electrically 
conducting fluid are considered. A scaling group of 
transformation is applied to the governing equations and the 
boundary conditions. 

The objective of present investigation is to study radiation 
effects on MHD flow over a porous vertical plate. To finding 
the solution, authors are using scaling group method of 
transformation. The attempt has also been made to study the 
effects of radiation, suction/ injection, porosity, thermal 
radiation, magnetic parameter, Schmidt number, Prandtl 
number on the fluid flow and the rate of heat and mass 
transfer. 


II. Mathematical Formulation of the Problem 


Consider a two dimensional steady laminar flow of an 
incompressible MHD fluid over a vertical porous plate in 
presence of constant suction and heat source and transverse 
magnetic field. Let the X-axis be taken in vertically upward 
direction along the plate and y -axis normal to it. The physical 
sketch and geometry of the problem is shown in Figure 1 . The 
fluid velocity u and V are along x and y -axis respectively. 


A magnetic field B is applied in the y -direction that is normal 
to the flow direction. Suction or injection is imposed on the 
porous plate. The temperature of the surface is held uniform at 
T w which is higher than the ambient temperature T.,. The 
species concentration at the surface is maintained uniform at 
C w = 1 while the ambient fluid concentration is assumed to be 
Co,. The governing fluid flow equations of continuity, 
momentum, heat transfer and mass transfer respectively are as 
follows: 
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porous vertical plate 


The boundary conditions are given by 
u = u w = x, v = — v w , C = C w , T = T w at y = 0 (5a) 

u — » 0,T — as y — > oo (5b) 

Here U and V are the velocity components in x and y 
directions, respectively, V is the kinematic coefficient of 
viscosity, [X is coefficient of viscosity, K is vortex viscosity, 
O is electrical conductivity of the fluid, K is porosity 
parameter, p is fluid density, T is fluid temperature, (3 T is 
thermal expansion coefficient, P c is concentration expansion 


coefficient, k is the thermal conductivity, C p is specific heat, 
D is mass diffusivity. 

By using Rosseland’s approximation, the radiative heat flux 
q r is given by 


4q, dT 4 

3kj 8y 


(6) 


Here Oj is the Stefan-Boltzman constant and k , is the 
absorption coefficient. It is assumed that the temperature 
variation within the flow is such that T 4 may be expanded in a 
Taylor series about T w and neglecting higher order terms, we 
get 

T 4 = 4T T 2 - 3T 4 
Equations ( 6 ) and (7) give 

8q r _ 160^ d 2 T 
8y 3kj dy 2 

Using equation ( 8 ) the energy equation (3) becomes 
_ ST _ ST _ k* S 2 T 16a<$* S 2 T p f SuY ( 9 ) 
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A. Nondimensionalization 

Introducing the following dimensionless variables 
considered by Hamad et al. (2012): 
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where Re = u w L/v is the Reynolds number, \|/ is the stream 
function, L being the characteristic length and u„o is reference 
velocity. 

Hence, equations (2), (9) and (4) reduce in the following 
form: 
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Hartmann number, K _ Up is the porosity parameter, 
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the Schmidt number. 


B. Group Transformation 

The application of group transformations has been 
considered to find similarity reduction of equations ( 11 ), ( 12 ) 
and (13). Consider the following group transformations 
r : x # = xn ai , y* = y tf 2 , y # = , (15) 

0 # = 0 EI a4 , tp # =(pil “ 5 

where ai, a 2 , a 3 , a 4 and a 5 are constants and Q is the 
parameter of point transformation. Now finding the relation 
among a’s such that 

A . (x # , y * , \|/ # , e # , cp # = Hj(x, y,\|/,0,(p, , • 5 ^ r ) 

dy* Sy 

^(x,y, V ,0,cp, ,^)(j = 1,2,3) 

By 3 

Ai, A 2 and A 3 are conformally invariant under the group 
transformation (15). 

By using above group transformation in equation (11) and 
solving the resulting equations one finds the following 
a 2 = a 4 =a 5 = 0 and = a 3 (16) 

Similarly equations (12), (13) and (14) are also giving 
a 1 =a 3 , a 2 = a 4 =a 5 =0, so these equations show 
invariant under the group transformation (15). 

Now the characteristic equations are 
dx _ dy _ dv|/ _ d 0 _ dtp ^ 7 ) 

x _ 0 _ \|/ _ 0 _ 0 


which give the following similarity transformations: 
r| = y, \|/ = xf(q), 0 = 0 (r|), and 9 = 9 ( 1 )) (18) 

Using these transformations, the momentum, energy and mass 
equations become 
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e » = z 3RPr( f9 , + Ecf „ 2 ) 

3R + 4 V ’ 

cp" = Sc(tpc — ftp') 

Subject to the boundary conditions 

f=f , f'=l, 0=1, cp = 1 at r| = 0 
w Y 1 


( 20 ) 

( 21 ) 


f'— >0, 9^0, cp— >0 as r| — > OO (22) 


The physical quantities of interest are the Skin friction 
coefficient Cf, Nusselt number Nu and Sherwood number Sh, 
which are defined as 


If), 




III. Method of solution 

The system of ordinary differential equations (19), (20) and 
(21) subject to the boundary conditions (22) have been solved 
numerically using Runge-Kutta method with shooting 
technique. The computations we are carried out using step size 
of Aq = 0.01 selected to be satisfactory for a convergence 
criterion of 10 6 in all cases. 

The physical quantities skin friction coefficient Cf, Nusselt 
number Nu and Sherwood number Sh indicate the wall shear 
stress, rate of heat transfer and rate of mass transfer 
respectively and these are proportional to the numerical values 
of f (0), -0 (0) and — <p'(0) respectively. 

IV. Results and discussion 

The numerical results for velocity, temperature, 
concentration, wall heat transfer, the rate of heat and mass 
transfer have been computed and represented. The 
computations are carried out for different values of different 
parameters for suction f w (= 0.5, 1), for injection f w (= -0.1, - 
0.5, -1). The calculations were made by taking radiation 
parameter (0.2 < R < 1), Sc = (0.22, 0.72, 6.8, 10) with Gr = 
0.1, Gc = 0.1, Ec = 1, Re = 1 and Prandtl number Pr = 0.72 for 
air at 1 atmospheric pressure. Different values of magnetic 
parameter, i.e. 0.1 < M < 1 were taken to analyse the study. 
The values for the skin friction coefficient, Nusselt number 
and Sherwood number have been tabulated in table 1. It is 
seen that the value of skin friction coefficient varies and is 
greater for suction while smaller for injection. Nusselt number 
is less for suction in comparison to injection and Sherwood 
number is higher for suction while smaller for injection. 


Table 1: Numerical values of f (0), 0 (0), rp (0) for different 
values of non dimensional parameters f w = 0.5, K = 0.1, Sc = 


R 

0.2 

0.5 


0.1 

0.1 

0.1 

0.4 

0.7 


0.5, Re = 1, Gr = 0.1, Gc = 0.1 and Ec = 1. 


M Pr 
0.5 0.72 

0.5 0.72 

0.5 0.72 

0.5 0.72 

0.5 0.72 


f "(0) 

0.097521 

0.047521 

0.000014 

0.352910 

0.838471 


0 ( 0 ) 

-0.149050 

-0.231590 

-0.295051 

-0.286591 

-0.241011 


<P'(0) 

-0.522674 

-0.515161 

-0.507571 

-0.530115 

-0.563746 


0.1 0.1 0.72 -0.452104 

0.1 1 0.72 0.659101 

0.1 0.5 0.22 0.075901 

0.1 0.5 6.8 0.010013 

0.1 0.5 10 0.001004 

0.1 0.5 0.72 0.000014 

0.1 0.5 0.72 0.001020 

0.1 0.5 0.72 0.001001 

0.1 0.5 0.72 0.015202 


-0.218399 -0.464637 

-0.334052 -0.551361 

-0.143201 -0.519201 

-1.327659 -0.511292 

-1.746591 -0.510820 

-0.286051 -0.504839 

-0.309251 -0.511870 

-0.298502 -0.689439 

-0.306302 -0.863251 


Fig. 2-4 exhibit the effects of suction on velocity, 
temperature and concentration. It is seen that the velocity, 
temperature and concentration decreases with an increase in 
suction parameter. Fig. 5-6 show the effect of injection on 
velocity and concentration profile. Velocity first increases 
then sharply decreases while concentration increases with 
rising as the value of injection. Fig.7-8 depict the effect of 
radiation parameter on velocity and temperature profile, it has 
been seen that velocity and temperature decreases with its 
increasing values (for R = 0.2, 0.5, 1). It is observed that 
velocity and temperature decrease while concentration 
increases with an increase in radiation parameter. Fig.9 shows 
the effect of magnetic parameter on velocity. Velocity 
increases as the value of M increases. The effect of Prandtl 
number Pr on temperature is represented by fig. 10. It is 
observed that the velocity sharply decreases with its increasing 
values (for Pr = 0.22, 6.8, 10). The effect of Schmidt number 
Sc on concentration is represented through fig. 11. It has been 
observed that concentration decreases as Schmidt number 
increases. 

x 
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Fig.7: Effect of radiation on velocity 




Fig. 8: Effect of suction on temperature 



Fig.5: Effect of injection on velocity 



Fig.9: Effect of magnetic parameter on velocity 



Fig. 10: Effect of Prandtl number on temperature 



Fig.6: Effect of injection on concentration 
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Fig. 1 1 : Effect of Schmidt number on concentration 
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